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Abstract. This paper shows in detail the application of a new stochastic approach for the characterization 
of surface height profiles, which is based on the theory of Markov processes. With this analysis we achieve 
a characterization of the scale dependent complexity of surface roughness by means of a Fokker-Planck 
or Langevin equation, providing the complete stochastic information of multiscale joint probabilities. The 
method is applied to several surfaces with different properties, for the purpose of showing the utility of 
this method in more detail. In particular we show evidence of the Markov properties, and we estimate the 
parameters of the Fokker-Planck equation by pure, parameter-free data analysis. The resulting Fokker- 
Planck equations are verified by numerical reconstruction of the conditional probability density functions. 
The results are compared with those from the analysis of multi-affine and extended multi-affine scaling 
properties which is often used for surface topographies. The different surface structures analysed here show 
in detail the advantages and disadvantages of these methods. 
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1 Introduction 

Among the great variety of complex and disordered sys- 
tems the complexity of surface roughness is attracting a 
great deal of scientific interest [1-5]. The physical and 
chemical properties of surfaces and interfaces are to a sig- 
nificant degree determined by their topographic structure. 
Thus a comprehensive characterization of their topogra- 
phy is of vital interest from a scientific point of view as 
well as for many applications [6-8]. 

Most popular methods used today for the character- 
ization of surface roughness are based on the concepts 
of self-affinity and multi-affinity, where the multifractal 
f{a) spectrum has been regarded as the most complete 
characterization of a surface [2,3,9,10]. One example of 
a measure of roughness which is commonly used in this 
context is the rms surface width w r (x) — ((h(x) — h) 2 )^ 2 , 
where h(x) is the measured height at point x, ( • ) r de- 
notes the average over an interval of length r around the 
selected point x, and h the mean value of h(x) in that 
interval. Thus the roughness is measured at a specific lo- 
cation x and over a specific scale r. Then a scaling regime 
of the ensemble average (w r ) in r, if it exists, is analyzed 
according to (w") ~ r^ a , usually a G Q. Here, (•) de- 
notes the mean over the available range in x. For a more 
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thorough introduction into scaling concepts we refer the 
reader to the literature, e.g. [2,3,9,10]. Terms concerning 
scaling concepts which are used in this paper are rapidly 
introduced in section 0] Here, we have to note the follow- 
ing points which concern stochastic aspects of roughness 
analysis: First, the ensemble average (w r ) must obey a 
scaling law as mentioned above, and second, the statis- 
tics of w r (x) are investigated over distinct length scales 
r, thus possible correlations between w r (x) and w r i(x) on 
different scales r, r' are not examined. 



In this paper we want to give a deeper introduction 
into a new approach to surface roughness analysis which 
has recently been introduced by us [11,12] and by [13]. 
This method is based on stochastic processes which should 
grasp the scale dependency of surface roughness in a most 
general way. No scaling feature is explicitly required, and 
especially the correlations between different scales r and 
r' are investigated. To this end we present a systematic 
procedure as to how the explicit form of a stochastic pro- 
cess for the r-evolution of a roughness measure similar to 
w r (x) can be extracted directly from the measured sur- 
face topography. This stochastic approach turns out to be 
a promising tool also for other systems with scale depen- 
dent complexity such as turbulence [11,14,15], financial 
data [16,17], and cosmic background radiation [18]. Also 
this stochastic approach has recently enabled the numer- 
ical reconstruction of surface topographies [13]. 
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Here we demonstrate our ansatz by analysing a num- 
ber of data sets from different surfaces. The purpose is 
to show extensively the utility of this method for a wide 
class of rough surfaces. The examples show different kinds 
of scaling properties which are, in addition, briefly anal- 
ysed. Among these examples is a collection of road sur- 
faces measured with a specially designed profile scanner. 
Preliminary results of the analysis of one of these surfaces 
have already been presented [12]. AFM measurement data 
from an evaporated gold film have already been analysed 
in an earlier stage of the method [11]. Since then, the 
method has been significantly refined and extended. Mea- 
surements of a steel crack surface were taken by confocal 
laser scanning microscopy (CLSM) [5]. 

As a measure of surface roughness we use the height 
increment [19] 

h r (x) :=h(x + r/2)-h(x-r/2) (1) 

depending on the length scale r. For other scale depen- 
dent roughness measures, see [3,10]. The height increment 
h r is used because its moments, which are well-known as 
structure functions (see section 0J), are closely connected 
with spatial correlation functions. Nevertheless, it should 
be pointed out that our method presented in the follow- 
ing could be easily generalized to any scale dependent 
measure, e.g. the above-mentioned w r (x) or wavelet func- 
tions [20]. As a new ansatz, h r is regarded as a stochastic 
variable in r. Without loss of generality we consider the 
process as being directed from larger to smaller scales. The 
focus of our method is the investigation of how the surface 
roughness is linked between different length scales. 

In the remainder of this paper we will first summarize 
in section El some central aspects of the theory of Markov 
processes which form the basis of our analysis procedure. 
Details concerning the measurement data are presented in 
sectional their scaling properties are analyzed in section^] 
The Markov properties of our examples are investigated 
in section [5] In section HJ we estimate for each data set 
the parameters of a Fokker-Planck equation. The ability 
of this equation to describe the statistics of h r in the scale 
variable r is then examined in section [7| followed by con- 
cluding remarks in section |H1 



2 Surface roughness as a Markov process 

Complete information about the stochastic process would 
be available from the knowledge of all possible n-point, or 
more precisely n-scale, joint probability density functions 
(pdf) p(h\, ri; h 2 , r 2 ; . . . ; h n , r n ) describing the probability 
of finding simultaneously the increments h\ on the scale 
r%, hi on the scale r 2l and so forth up to h n on the scale 
r n . Here we use the notation hi(x) = h n (x), see eq. 
Without loss of generality we take r\ < r 2 < ... < r n . 
As a first question one has to ask for a suitable simplifi- 
cation. In any case the n-scale joint pdf can be expressed 
by multiconditional pdf 

p(hx,ri;...;h n ,r n ) = 



p(hi,rx\h 2 ,r 2 ; ■ ■ ■;h n ,r n ) ■ p(h 2 ,r 2 \h 3 ,r 3 ; . . .;h n ,r n ) 
■ ... ■ p(h n -i,r n -i\h n ,r n ) -p(h n ,r n ) . (2) 

Here, p(hi,ri\hj,rj) denotes a conditional probability of 
finding the increment hi on the scale under the condi- 
tion that simultaneously, i.e. at the same location x, on a 
larger scale rj the value hj was found. It is defined with 
the help of the joint probability p(hi, r*; hj, rj) by 

P( h *> r *\ h i> r >) = -q—j . (3) 
An important simplification arises if 
p{h u n\h i+ t,ri + i; ...;h n , r„) = p{h t , ri\h i+ i,r l+1 ) . (4) 

This property is the defining feature of a Markov process 
evolving from r i+1 to r^. Thus for a Markov process the 
n-scale joint pdf factorize into n conditional pdf 

p(hi,ri; h n ,r n ) =p{h\,ri \h 2 ,r 2 ) 

■ ... -p(h n -i,r n -i \h n ,r n ) -p(h n ,r n ) . (5) 

The Markov property implies that the r-dependence of 
h r can be regarded as a stochastic process evolving in r, 
driven by deterministic and random forces. Here it should 
be noted that if condition |TJJ holds this is true for a pro- 
cess evolving in r from large down to small scales as well 
as the reverse from small to large scales [21]. Equation 
also emphasizes the fundamental meaning of conditional 
probabilities for Markov processes since they determine 
any n-scale joint pdf and thus the complete statistics of 
the process. 

For any Markov process a Kramers-Moyal expansion of 
the governing master equation exists [22]. For our height 
profiles it takes the form 

d 

-r — p{h r7 r\h ,r ) = 

OO / \ k 

j^\dhr) Dik} ( h r^)P( h r^\ho,r a ). (6) 

The minus sign on the left side of eq. © expresses the 
direction of the process from larger to smaller scales, fur- 
thermore the factor r corresponds to a logarithmic vari- 
able p = In r which leads to simpler results in the case of 
the scaling behaviour [23]. To derive the Kramers-Moyal 
coefficients D^ k ' (h r ,r), the limit Ar — > of the condi- 
tional moments has to be performed: 

D (k) {h r ,r) = lim M (k) {h r ,r,Ar) , where (7) 

Ar^O 

M {k) (h r ,r, Ar) = 
-^,J + _ (h-K) k p(kr-Ar\h r ,r)dh. (8) 

The moments M^ k \h r ,r, Ar) characterize the alteration 
of the conditional probability p{h r , r\ho, ro) over a finite 
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step size Ar = r$ — r and are thus also called "transitional 
moments" . 

A second major simplification is valid if the noise in- 
cluded in the process is Gaussian distributed. In this case 
the coefficient Z)' 4 ' vanishes (from eqs. Q) and JSJ) it 
can be seen that D^' is a measure of non-gaussianity of 
the included noise). According to Pawula's theorem, to- 
gether with Z?( 4 ) all the D^ k > with k > 3 disappear and 
the Kramers-Moyal expansion JJ]J collapses to a Fokker- 
Planck equation [22], also known as Kolmogorov equa- 
tion [24]: 

d 

p{h r ,r\ho,r Q ) = (9) 

{~lk D{1){hr ' r) + ^2 D{2) (hr,r)}p(h r ,r\ho,ro) 

The Fokker-Planck equation then describes the evolution 
of the conditional probability density function from larger 
to smaller length scales and thus also the complete n-scale 
statistics. The term D^> (h r , r) is commonly denoted as 
the drift term, describing the deterministic part of the pro- 
cess, while D^ 2 \h r ,r) is designated as the diffusion term, 
determined by the variance of a Gaussian, <5-correlated 
noise (compare also eqs. and (JHJ). 

By integrating over ho it can be seen that the Fokker- 
Planck equation Q is also valid for the unconditional 
probabilities p(h r ,r) (see also section 0. Thus it covers 
also the behaviour of the moments (also called struc- 
ture functions) including any possible scaling behaviour. 
An equation for the moments can be obtained by addi- 
tionally multiplying with ft™ and integrating over ft r 

-r^{K) = (10) 
or 

n{D^\h r ,r)h^- x ) +n(n~l)(D^(h r ,r)h r r 1 - 2 ) . 

For D^> being purely linear in h r (D^ 1 ' = ah r ) and 
purely quadratic (D^ = Ph 2 .), the multifractal scaling 
(ft™) ~ r^" with £„ = na + n(n — l)/3 is obtained from 
H10(l . If in contrast is constant in h r , a monofractal 
scaling where are linear in n may occur, see [11]. 

Lastly, we want to point out that the Fokker-Planck 
equation (0 corresponds to the following Langevin equa- 
tion (we use Ito's definition) [22] 

- = D^(h r ,r)/r+ y/DW(h r ,r)/rr(r) , (11) 

where r(r) is a Gaussian distributed, (5-correlated noise. 
The use of this Langevin model in the scale variable opens 
the possibility to directly simulate surface profiles with 
given stochastic properties, similar to [13]. 

With this brief summary of the features of stochastic 
processes we have fixed the scheme from which we will 
present our analysis of diverse rough surfaces. There are 
three steps: First, the verification of the Markov property. 
Second, the estimation of the drift and diffusion coeffi- 
cients and D^ 2 \ Third, the verification of the es- 
timated coefficients by a numerical solution of the corre- 
sponding Fokker-Planck equation, thus reconstructing the 
pdf which are compared to the empirical ones. 



3 Measurement data 

With the method outlined in section [3 we analysed a col- 
lection of road surfaces measured with a specially designed 
profile scanner as well as two microscopic surfaces, namely 
an evaporated gold film and a crack surface of a low- 
alloyed steel sample, as already mentioned in section ^ 

The road surfaces have been measured with a specially 
designed surface profile scanner. The Longitudinal resolu- 
tion was 1.04 mm, the profile length being typically 20 m 
or 19000 samples, respectively. Between ten and twenty 
parallel profiles with a lateral distance of 10 mm were 
taken for each surface, see fig. ^ The vertical error was 
always smaller than 0.5 mm but in most cases approxi- 
mately 0.1mm. Details can be found in [25]. 

For the Au film data, the surface of four optical glass 
plates had been coated with an Au layer of 60 nm thick- 
ness by thermal evaporation [11]. The topography of these 
films was measured by atomic force microscopy at differ- 
ent resolutions, resulting in a set of images of 256 x 256 
pixels each, where every pixel specifies the surface height 
relative to a reference plane, see fig.|U Out of these images 
99 could be used for the analysis presented here, resulting 
in about 6.5 • 10 6 data points. Sidelengths vary between 
36 nm and 2.8 /im. 

The sample of the crack measurements was a fracture 
surface of a low-alloyed steel (german brand 10MnMoNi5- 
5). A detailed description of the measurements can be 
found in [5]. Three CLSM (Confocal Laser Scanning Mi- 
croscopy) images of size 512x512 pixels in different spa- 
tial resolutions were available, see fig. El Pixel sizes are 
0.49, 0.98, and 1.95 /im, resulting in image widths of 251, 
502, and 998 fim, repectively. Unavoidable artefacts of the 
CLSM method were removed by simply omitting for each 
image those data with the smallest and largest height 
value, similar to [5]. Nevertheless this cannot guarantee 
the detection of all the artefacts. The possible consequences 
are discussed together with the results. 

For the analysis in the framework of the theory of 
Markov processes, we will normalize the measurement data 
by the quantity defined by 

al = lim (h 2 r ) . (12) 

r — >oo 

Thus it is possible to obtain dimensionless data with a 
normalization independent of the scale r, in contrast to 
e.g. a 2 = (hi). As a consequence the results, especially 
M^ k \h ri r,Ar) and D^ k \h r ,r) (cf. section©, will also be 
dimensionless. From the definition it is easy to see that 
(Too can be derived via a 2 ^ = 2o 2 x = 2((h(x) — h) 2 ) if h(x) 
becomes uncorrelated for large distances r. 



4 Scaling analysis 

In this paper a number of examples was selected from 
all data sets under investigation. Because most popular 
methods of surface analysis are based on scaling features 
of some topographical measure, the examples were chosen 
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with respect to their different scaling properties as well 
as their results from our analysis based on the theory of 
Markov processes. 

In the analysis presented here we use the well-known 
height increment h r (x), which has been defined in eq. ^ 
as a scale-dependent measure of the complexity of rough 
surfaces [19]. Scaling properties are reflected by the r- 
dependence of the so-called structure functions 

S n (r) = (\h?\) . (13) 

If one then finds 

S n (r) ~ r*» (14) 

for a range of r, this regime is called the scaling range. In 
that range the investigated profiles have self-affine proper- 
ties, i.e., they are statistically invariant under an anisotropic 
scale transformation. If furthermore the dependence of the 
exponents on the order n is nonlinear, one speaks of 
multi-affine scaling. Those properties are no longer iden- 
tified by a single scaling exponent, but an infinite set of 
exponents. A detailed explanation of self- and multi-affine 
concepts is beyond the scope of this article. Instead, we 
would like to refer the reader to the literature [2,3,9, 10]. 
The power spectrum, which often is used to determine 
scaling properties, can easily be derived from the second 
order structure function. It is defined as the Fourier trans- 
form of the autocorrelation function R(r), which itself is 
closely related to S 2 {r) by R(r) = (h(x) 2 ) - S 2 (r)/2, by 
comparing eqs. JD and (fH^ . 

In addition to the r-dependence of the structure func- 
tions, a generalized form of scaling behaviour can be deter- 
mined analogously to the Extended Self Similarity (ESS) 
method which is popular in turbulence research [26] . When 
the S n (r) are plotted against a structure function of spe- 
cific order, say S 3 (r), in many cases an extended scaling 
regime is found according to 

S n (r)~(S 3 (r)) U . (15) 

Clearly, meaningful results are restricted to the regime 
where S 3 is monotonous. It is easy to see that now the £ n 
can be obtained by 

£n = Cn ' & , (16) 

cf. [26] . While for turbulence it is widely accepted that by 
this means experimental defiencies can be compensated 
to some degree, for surface roughness the meaning of ESS 
lies merely in a generalized form of scaling properties. 

It should be noted that the results of any scaling anal- 
ysis may be influenced by the method of measurement, 
by the definition of the roughness measure, here h r (x) (or 
w r (x) as mentioned in section as well as by the al- 
gorithms used for the analysis [5,27]. Nevertheless, this 
problem is not addressed here as the main focus of our 
investigations is the application of the theory of Markov 
processes to experimental data. 



4.1 Surfaces with scaling properties 

In fig. ^ we present road surface data with different kinds 
of scaling properties. For each data set a short profile sec- 
tion is shown. Structure functions of order one to six on 
double logarithmic scale are presented in fig. [21 Follow- 
ing the arguments in [14], higher order structure func- 
tions cannot be evaluated with sufficient precision from 
the given amount of data points. The worn asphalt pave- 
ment (Road 1) is an example of a comparably large scaling 
regime over more than one order of magnitude in r. A sur- 
face with similar features, namely a cobblestone road, has 
already been presented in [12]. Two separate scaling re- 
gions are found for a Y-shaped concrete stone pavement 
(Road 2). Additionally a sharp notch can be seen in the 
structure functions at r = 0.2 m, indicating a strong peri- 
odicity of the pavement caused by the length of the indi- 
vidual stones. The third example, a "pebbled concrete" 
pavement (Road 3), consists of concrete stones with a 
top layer of washed pebbles. This material is also known 
as "exposed aggregate concrete" . Here, the scaling region 
of the structure functions is only small. For the basalt 
stone pavement (Road 4), being the fourth example, scal- 
ing properties are poor. We have nevertheless marked a 
possible scaling range and derived the respective scaling 
exponents for comparison with the other examples. Simi- 
lar to the Y-shaped concrete stones, a periodicity can be 
found at about 0.1 m length scale. 

The results for the generalized scaling behaviour ac- 
cording to eq. I|15|) are shown in fig. O It can be seen that 
indeed for three of the surfaces in fig. ^ an improved scal- 
ing behaviour is found by this method. Only for Road 4 do 
the generalized scaling properties remain weak. In fig. 0] 
the scaling exponents £ n of the structure functions within 
the marked scaling regimes in fig. [21 were determined and 
plotted against the order n as open symbols. Additionally, 
values of £„ were derived according to eqs. \ib\ and (|16fl 
and added as crosses. For Road 2 two sets of exponents 
correspond to the two distinct scaling regimes in fig. [3 
Even though there is only one set of ( n found in fig. [31 
two sets of are obtained due to the two different £3, see 
eqs. JUJl and l(Tf])l. 

All surfaces from fig. ^ show a more or less nonlinear 
dependence of the on n, indicating multi-affine scal- 
ing properties. The scaling exponents obtained via the 
generalized scaling according to eq. (|TS)l are in good cor- 
respondence with the achieved by the application of 
eq. (|14|l . Deviations are seen for Road 2 and at higher 
orders for Road 3, possibly caused by inaccuracies in the 
fitting procedure. For Road 4 no generalized scaling is ob- 
served (compare fig. [^ and thus values of £ n cannot be 
derived from £„. From this we conclude that scaling prop- 
erties for some cases are questionable as a comprehensive 
tool to characterize the complexity of a rough surface. 

An example for good scaling properties is the gold film 
surface (Au). To increase statistical accuracy, increments 
are evaluated here in the direction of the rows of the im- 
ages as well as the columns. In fig. [3 two of the 99 images 
under investigation are shown. Figure[S|presents the struc- 
ture functions S n (r), derived from all images. The surface 
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Fig. 1. Measurement data from selected road surfaces with 
different kinds of scaling properties. Pavements are worn as- 
phalt (Road f), Y-shaped concrete stones (Road 2), pebbled 
concrete stones (Road 3), and basalt stones (Road 4), from top 
to bottom. For each surface a short section of the respective 
height profile is shown. 



Fig. 2. Structure functions S n (r) of selected road surfaces 
(see fig. with different kinds of scaling properties on a log- 
log scale (see text). Symbols correspond to orders n; diamonds 
(n = 1), triangles (n = 2), circles (n = 3), squares (n = 4), x 
signs (n = 5), and plus signs (n = 6). 
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Fig. 3. Generalized scaling analysis of the surfaces shown in Fig. 4. Scaling exponents £ n of the surfaces shown in fig. 
fig. Q Structure functions S n are displayed versus S 3 on a log- achieved via eq. 1141 (open symbols) and those obtained via 
log scale. Symbols correspond to orders n as in fig. [5] £ n from eq. 116H (crosses) . For Road 2 two sets of exponents 

were obtained from the two scaling regimes found for S n (r) in 

fig- HI 
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is randomly covered with granules which show no typical 
diameter. A scaling regime of more than one order of mag- 
nitude in r is found for the structure functions S n (r) in 
fig. El Generalized scaling behaviour is clearly present as 
shown in fig-E^a). The scaling exponents £ n presented in 
part (b) of the same figure are nearly linear in n, thus this 
surface can not be regarded as multi-affine, but appears 
to be self-affine. Here, the achieved via eq. (|16fl match 
perfectly those obtained from eq. 1|15|) . 





50 nm 0.5 fim 

Fig. 5. AFM images of the Au film surface. Sidelengths are 
110 nm and 1.1 fjxci. The relative surface height is represented 
as gray level. Maximum heights are 7.2 nm and 13.3 nm, re- 
spectively. 
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Fig. 7. Generalized scaling properties (a) and scaling expo- 
nents (b) of the Au film surface shown in fig. Scaling expo- 
nents £ n achieved via eq. 115H are marked by open symbols, 
those obtained via £ n from eq. 1 Kit by crosses. Compare also 
with fig. 
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Fig. 6. Structure functions S n (r) of the Au film surface on a 
log-log scale. Symbols correspond to orders n as in fig. El 



4.2 Surfaces without scaling properties 

To complete the set of examples, we present two surfaces 
without scaling properties. The first one is a smooth as- 
phalt road (Road 5), shown in fig. |Hl No power law can 
be detected for the S n (r) but a generalized scaling is ob- 
served in fig. lTTT ab The range of values of S* 3 , however, is 
relatively small. 

The second example lacking a scaling regime is the 
steel fracture surface (Crack). One of the three CLSM im- 
ages under investigation is shown in fig-Efa). Figure EJb) 
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Fig. 8. Measurement data (a) and structure functions (b) from 
a road surface without scaling properties (Road 5). The pave- 
ment is smooth asphalt. See also figs. and |5] 
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presents an additional REM image at a higher resolution, 
which gives an impression of the surface morphology. For 
the structure functions in fig. 1101 no scaling properties are 
found, and the dependences of S n (r) on 5 3 (r) in ng. lTlT b) 
also deviate from proper power laws. It should be noted 
that in general scaling properties not only depend on the 
respective data set but also on the analysis procedure. 
Using other measures than h r (x), in [5] scaling regimes 
of those measures have been found, and scaling exponents 
could be obtained. 



(a) 



(b) 























' 200 /iin 


' 100 //m 1 



Fig. 9. Measurement data from a steel crack surface (Crack), 
(a) CLSM image, side length 502 /im, (b) REM image, side 
length 140 /jm. 
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Fig. 11. Generalized scaling properties of the surfaces shown 
in figs. |H1 (Road 5 (a)) and (Crack (b)). Compare also with 
figs. H and □ 
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Fig. 10. Structure functions S n (r) of the CLSM images from 
a steel crack surface (Crack) on a log-log scale. The symbols 
correspond to orders n as in fig.|21 



As outlined in sections ^ an d 12 we want to describe the 
evolution of the height increments h r {x) in the scale vari- 
able r as realizations of a Markov process with the help 
of a Fokker-Planck equation. Consequently, the first step 
in the analysis procedure has to be the verification of the 
Markov properties of h r (x) as a stochastic variable in r. 

For a Markov process the denning feature is that the 
n-scale conditional probability distributions are equal to 
the single conditional probabilities, according to eq. QJ. 
With the given amount of data points the verification of 
this condition is only possible for three different scales. 
Additionally the scales r are limited by the available pro- 
file length. For the sake of simplicity we will always take 
r% — r 2 = r 2 — r% = Ar. Thus we can test the validity of 
eq. (0J in the form 



4.3 Conclusions on scaling analysis 

To conclude the scaling analysis of our examples, we have 
chosen surfaces with a range of different scaling properties 
from good scaling over comparably wide ranges, such as 
for Road 1 and Au, to the absence of scaling, such as 
for Road 5 and Crack. The generalized scaling analysis, 
analogous to ESS [26] , leads to the same scaling exponents 
as the dependence of the structure functions S n (r) on the 
scale r, with some minor deviations. 



p{h x , n | h 2 , n+Ar) = p(h 1 ,r 1 \h 2 , n+Ar; h 3 = 0, n+2Ar) . 

(17) 

Note that in eq. (|17fl we take /13 — to restrict the number 
of free parameters in the pdf with double conditions. 

Three procedures were applied to find out if Markov 
properties exist for our data. From the results of all three 
tests we will find a minimal length scale Im for which 
this is the case. The meaning of this so-called Markov 
length will be discussed below. In the following we will 
demonstrate the methods using the example of the Au 
surface. 
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5.1 Testing procedures 

The most straightforward way to verify eq. I|17f) is the vi- 
sual comparison of both sides, i.e., the pdf with single and 
double conditions. This is illustrated in fig.Elf 01 ' two dif- 
ferent scale separations Ar = 17nm and 35 nm. In each 
case a contour plot of single and double conditional proba- 
bilities p{h 1 ,r 1 \h 2 ,r 2 ) and p(hi, r\\h 2 , r 2 ; /&3=0, r 3 ) is pre- 
sented in the top panel of (a) and (b), respectively. Below 
two one-dimensional cuts at fixed values of h 2 « icr^ 
are shown, representing directly p(h%, ri\h 2 =±a ca , r 2 ; h 3 — 
0, ra). It can be seen that in panel (a), for the smaller value 
of Ar, the single and double conditional probability are 
different. This becomes clear from the crossing solid and 
broken contour lines of the contour plot as well as from the 
differing lines and symbols of the one-dimensional plots 
below. Panel (b), for Ar = 35 nm, shows good correspon- 
dence of both conditional pdf. We take this finding as a 
strong hint that for this scale separation Ar eq. I@J is valid 
and Markov properties exist. Following this procedure for 
all accessible values of Ar, the presence of Markov prop- 
erties was examined. For this surface Markov properties 
were found for scale distances from (25 ± 5) nm upwards. 

The validity of eq. I)17|l can also be be quantified math- 
ematically using statistical tests. An approach via the 
well-known x 2 measure has been presented in [28] , whereas 
in [14] the Wilcoxon test has been used. Next, we give a 
brief introduction to this procedure, which will be used 
here, too. More detailed discussions of this test can be 
found in [14,21,29]. For this procedure, we introduce the 
notation of two stochastic variables xi,i — 1, . . . ,n and 
Uj , j = 1 , . . . , to which represent the two samples from 
which both conditional pdf of eq. are estimated, i.e. 



x(h r2 ,n,r 2 ) = h ri \hr 2 
y(hr 2 ,h r3 ,ri,r 2 ,r 3 ) = h ri | h . hr 



(18) 



Here denotes the conditioning. All events of both sam- 
ples are sorted together in ascending order into one se- 
quence, according to their value. Now the total number of 
so called inversions is counted, where the number of in- 
versions for a single event yj is just the number of events 
of the other sample which have a smaller value xt < yj . If 
eq. (|17fl holds and n,m > 25, the total number of inver- 
sions Q is Gaussian distributed with 



(Q) = nm/2 and 
<7q = \J nm(n + m + 1) /12 



(19) 



We normalize Q with respect to its standard deviation 
and consider the absolute value 



t = \Q-(Q)\/<T Q 



(20) 



For its expectation value it is easy to show that (t) = 
\TlpK (still provided that l|17|) is valid), where here the 
average (•) is performed over h 2 . If a larger value of (t) is 
measured for a specific combination of r and Ar, we con- 
clude that eq. l(T7J) is not fulfilled and thus Markov proper- 
ties do not exist. A practical problem with the Wilcoxon 
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Fig. 12. Test for Markov properties of Au film data for two 
different scale separations Ar = 14 nm (a) and 35 nm (b), 
where Ar — rz — r-z = r-z — n (see text). In both cases r2 = 
169 nm. In each case a contour plot of conditional probabili- 
ties p(hi, T\ |/i2, 5*2) (dashed lines) and p(hi, ri|/i2, r2\ /i3=0, rg) 
(solid lines) is shown in the top panel. Contour levels differ 
by a factor of 10, with an additional level at p = 0.3. Be- 
low the top panels in each case, two one-dimensional cuts at 
/12 ~ ±o"oo are shown with p(hi, ri|/i2, ^2) as dashed lines and 
p(h 1 ,r 1 \h2,r2; h3=0,rz) as circles. 
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test is that all events Xi,yj have to be statistically in- 
dependent. This means that the intervals of subsequent 
height increments h r have to be separated by the largest 
scale involved. Thus the number of available data is dra- 
matically reduced. 



5 p 

4 - 

3 - 
2 

1 : 




(a) 



10 20 30 40 50 60 70 
Ar [nm] 

Fig. 13. Wilcoxon test for the Au surface. The scale r is 28 nm. 
The theroretically expected value (t) — y2/7r is marked with 
a horizontal line, the Markov length Im = (25 ± 5) nm with a 
vertical line. 



In fig. H3l we present for the Au surface measured values 
of (t(r, Ar)) at a scale r = 28 nm. The Markov length Im 
is marked where (t(r, Ar)) has approached its theroretical 
value a/2/-7t. 

Another method to show the validity of condition 117fl 
is the investigation of the well-known necessary condition 
for a Markovian process, the validity of the Chapman- 
Kolmogorov equation [22] 



p(hi,ri\h 3 ,r 3 ) = 



+ OO 



p(hi,ri\h,2, r 2 )p(h 2 , r 2 \h 3 , r 3 )dh 2 ■ 

(21) 

We use this equation as a method to investigate the Markov 
properties of our data. This procedure was used for exam- 
ple in [13,30-32] for the verification of Markov properties. 
It also served to show for the first time the existence of 
a Markov length in [28]. The conditional probabilities in 
eq. (|21|l are directly estimated from the measured profiles. 
In fig. El both sides of the Chapman-Kolmogorov equa- 
tion are compared for two different values of Ar. In an 
analogous way to fig. ^1 for each Ar the two conditional 
probabilities are presented together in a contour plot as 
well as in two cuts at fixed values of h 2 . While for the 
smaller value Ar = 14 nm both the contour lines and the 
cuts at fixed h 2 clearly differ, we find a good correspon- 
dence for the larger value Ar = 35 nm. 

A third method which we did not use here but which 
is reported in the literature is based on the description of 
the stochastic process by a Langevin equation. With this 
knowledge of the Langevin equation l|ll|) the noise can be 
reconstructed and analyzed with respect to its correlation 
[33,34]. 




3-2-1 12 3 
hi [a J 
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Fig. 14. A check of the Chapman-Kolmogorov equation l|21|l 
for the Au surface for two different scale separations Ar = 
14 nm (a) and 35 nm (b). In both cases r2 = 169 nm. The plots 
are organized in the same way as fig. 1121 The pdf representing 
the left side of J2TJ 

are shown with solid lines, the integrated 
pdf of the right side of 1211 as dashed lines and circles. 

5.2 Conclusions on Markov properties 

The results of the methods described above were combined 
to determine whether Markov properties of the height in- 
crement h r {x) in the scale variable are present for our 
surface measurements. We found Markov properties for 
all the selected examples of surface measurement data. It 
is also common to all examples that these Markov prop- 
erties are not universal for all scale separations Ar but 
there exists a lower threshold which we call the Markov 
length Im- It was determined in each case by systematic 
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Surface 


Im 


Surface 


Im 


Surface 


Im 


Road 1 


33 mm 


Road 2 


10 mm 


Road 3 


4.2 mm 


Road 4 


17 mm 


Road 5 


4.2 mm 


Au 


25 nm 


Crack 


20 fim 











Table 1. Markov lengths Im for all the surfaces presented. 



application of the three testing procedures for all accessi- 
ble length scales r and scale separations Ar. The resulting 
values are listed in table The presence of Markov prop- 
erties only for values of Ar above a certain threshold has 
also been found for stochastic data generated by a large 
variety of processes and especially occurs in turbulent ve- 
locities [11,14,18,28,35]. 

The meaning of this Markov length I m may be seen in 
comparison with a mean free path length of a Brownian 
motion. Only above this mean free path is a stochastic 
process description valid. For smaller scales there must be 
some coherence which prohibits a description of the struc- 
ture by a Markov process. If for example the description 
of a surface structure requires a second order derivative in 
space, a Langevin equation description l|ll|) becomes im- 
possible. In this case a higher dimensional Langevin equa- 
tion (at least two variables) is needed. It may be interest- 
ing to note that the Markov length we found for the Au 
surface of about 30 nm coincides quite well with the size 
of the largest grain structures we see in fig. [S] Thus the 
Au surface may be thought of as a composition of grains 
(coherent structures) by a stochastic Markov process. 

In the case of Road 2 with its strong periodicity at 
0.2 m the Markov properties end slightly above this length 
scale. It seems evident that here the Markov property is 
destroyed by the periodicity. While some of the other sur- 
faces also have periodicities, these are never as sharp as 
for Road 2. An upper limit for Markov properties could 
not be found for any of the other surfaces. 

Another interesting finding can be seen from figs. P21 
El El an d El There is no connection between the scal- 
ing range and the range where Markov properties hold. 
Regimes of scaling and Markov properties are found to be 
distinct, overlapping or covering, depending on the sur- 
face. Data sets which fulfill the Markov property do not 
in all cases show a scaling regime at all. Also, on the other 
hand, scaling features seem not to imply Markov proper- 
ties, which has been indicated previously for some numer- 
ically generated data in [11]. While there is always an up- 
per limit of the scaling regime, we found only one surface 
for which the Markov properties possess an upper limit. 

6 Estimation of drift and diffusion coefficients 

As a next step we want to concentrate on extracting the 
concrete form of the stochastic process, if the Markov 
properties are fulfilled. As mentioned in sectionFJlour anal- 
ysis is based on the estimation of Kramers-Moyal coeffi- 
cients. The procedure we use to obtain the drift (D^ 1 ') and 
diffusion coefficient (D^) for the Fokker-Planck equa- 
tion © was already outlined by Kolmogorov [24] , see also 



[14,22]. First, the conditional moments M^ k \h r , r, Ar) for 
finite step sizes Ar are estimated from the data via the 
moments of the conditional probabilities. This is done by 
application of the definition in eq. JBJ), which is recalled 
here: 

M^(h r ,r,Ar) = J + {h-h r ) k p(h,r-Ar\h r ,r) dh 

(22) 

The conditional probabilities in the integral are obtained 
by counting events in the measurement data as shown al- 
ready in sectional Here, one fundamental difficulty of the 
method arises: For reliable estimates of conditional prob- 
abilities we need a sufficient number of events even for 
rare combinations of h, h r . Consequently, a large amount 
of data points is needed. This problem becomes even more 
important if one takes into account that a large range in 
r should be considered. The number of statistically inde- 
pendent intervals h r is limited by the length of the given 
data set and decreases with increasing r. 

In a second step, the coefficients D^ k \h r ,r) are ob- 
tained from the limit of M^ k ' {h r , r, A r ) when Ar approaches 
zero (see definition in eq. 0). For fixed values of r and h r 
a straight line is fitted to the sequence of (h r ,r, Ar) 
depending on Ar and extrapolated against Ar = 0. The 
linear dependence corresponds to the lowest order term 
when the Z\r-dependence of M^ k \h r ,r, Ar) is expanded 
into a Taylor series for a given Fokker-Planck equation 
[33,36]. Our interpretation is that this way of estimating 
the is the most advanced one, and also performs bet- 
ter than first parameterizing the M k and then estimating 
the limit Ar — > for this parameterization, as previously 
suggested in [14,21]. 

There have been suggestions to fit functions to M^ k ' 
other than a straight line, especially for the estimation of 
D^ 2 \ see [21]. Furthermore it has been proposed to use 
particular terms of the above-mentioned expansion to di- 
rectly estimate D^ k \h r ,r) without extrapolation [32]. On 
the other hand, in [37] it becomes clear that there can be 
manifold dependences of on Ar which in general are 
not known for a measured data set. Consequently, one may 
state that there is still a demand to improve the estimation 
of D^ k \ At the present time we suggest to show the qual- 
ity of the estimated by verification of the resulting 
Fokker-Planck equation, once its drift and diffusion coeffi- 
cients have been estimated. However, for our data neither 
nonlinear fitting functions nor correction terms applied to 
the resulted in improvements of the estimated . 

A crucial point in our estimation procedure is the range 
of Ar where the fit can be performed. Only those Ar can 
be used where Markov properties were found in the scale 
domain. In section|3]we showed that for our data Markov 
properties are given for Ar larger than the Markov length 
Im (see table ^ • In order to reduce uncertainty, a large 
range of Ar as the basis of the extrapolation is desirable. 
From eq. iTJ^J it can be seen, however, that Ar must be 
smaller than r. As a compromise between accuracy and 
extending the scale r to smaller values, in many cases an 
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extrapolation range of Im < Ar < 21m was used (cf. table 
|2J) • This procedure is shown in fig. El for Road 1. 



Surface 


hi 


Imax 


Surface 


Im 


Imax 


Road 1 


33 mm 


67 mm 


Road 2 


10 mm 


21 mm 


Road 3 


4.2 mm 


19 mm 


Road 4 


17 mm 


25 mm 


Road 5 


4.2 mm 


8.3 mm 


Au 


25 nm 


84 nm 


Crack 


20 /im 


44 jam 









Table 2. Extrapolation ranges for all the presented surfaces. 
Listed are the smallest (Im) and largest (l ma x) values of Ar 
used for extrapolation of the D^ k \h r ,r). 
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Fig. 15. Extrapolation procedure for (a) and (b), 

illustrated for surface Road 1. Length scale r is 108 mm, h r 
is — <Joo- Values of M' 1 ' and M' 2 ' inside the range marked 
by broken lines were used for the extrapolation. The results 
D w (h r ,r) and D (2) (/i r ,r) are marked with filled circles on 
the vertical axis. 



6.1 Estimation results 

Following the procedure outlined above, D^(h r ,r) and 
D^ 2 \h r ,r) were derived for the measurement data pre- 
sented in section [31 with the exception of Road 5 (see 
below). For the road surfaces, estimations were performed 
for length scales r separated by ten measurement steps 
or 10.4 mm, respectively, to reach a sufficient density over 
the range where the coefficients were accessible. FigureslTBI 



and 1 171 show estimations of the drift coefficients and 
the diffusion coefficients D 1 - 2 ^ for the road surfaces, each 
performed for one fixed length scale r. The error bars are 
estimated from the errors ofMW(/i r ,r, Ar) via the num- 
ber of statistically independent events contributing to each 
value, assuming that each bin ofp(hi, fil^o, r o) containing 
N events has an intrinsic uncertainty of ±^/N. Addition- 
ally, values of are added to the plots of which 
have been estimated in the same way. Thus it can be seen 
that in all cases is small compared to D^ 2 \ except for 
Road 4, and in most cases its statistical errors are larger 
than the values themselves. Negative values are not shown 
because the vertical axes start at zero. As is posi- 

tive by definition, the occurence of negative values of 
results from the limit Ar — > and should be only due to 
the statistical errors involved. Even if there is no evidence 
that Z?( 4 ) is identically zero, the presented values give a 
hint that its influence in the Kramcrs-Moyal expansion 
© is rather small and the assumption of a Fokker-Planck 
equation 10 is justifiable, with the possible exception of 
Road 4. 

Estimated drift and diffusion coefficients and 
for the Au surface are shown in fig. ^| for r — 169 nm. 
Again, was added to the plot of D^ 2 \ in this case 
without error bars to enhance clarity. Errors of are 
in this case always much larger than the values themselves 
and would cover the values of D^> as well as their er- 
rors. Also the error bars of D^> appear to be quite large 
for the Au surface. The data here are measured as two- 
dimensional images, thus the number of statistically inde- 
pendent h r (x) decreases quadratically with increasing r, 
resulting in rather large error estimates. For the calcula- 
tion of Z)( fc ) nevertheless all accessible h r (x) were used. 
As the regime of Markov properties starts at Ar = 25 nm, 
the range 25 nm < Ar < 84 nm was used as basis for 
the extrapolation (see table |2J- For r < 84 nm the upper 
limit was reduced in order to derive the coefficients also 
for smaller scales r (compare also section |HJ). In this way 
the drift and diffusion coefficients of the Au film could be 
worked out from 281 down to 56 nm. 

In the same way as for the other surfaces, estimations 
of the Kramers-Moyal coefficients were performed for the 
steel crack. The results are shown in fig. ^3 Again, the 
estimates for are also presented, which are of the 

same order of magnitude as for higher values of h r 
(\h r \ > 0.5croo). 

For the surface Road 5 (cf. fig. drift and diffusion 
coefficients could not be estimated. The reason can be 
seen in fig. [2U] The diagram shows the dependence of 
M^(h r ,r,A r ) and M (2 \h r ,r, Ar) on Ar for fixed r and 
h r , in this case 104mm and —0.6(700. For Ar > Im it 
can be seen that and M (2) behave like 1/Ar. This 

behaviour can be explained by the presence of some ad- 
ditional uncorrelated noise, where additional means inde- 
pendent of the stochastic process. A similar behaviour was 
found for financial market data [38]. In this case the in- 
tegral in eq. (|22J) will tend to a constant for small Ar, 
independent of the value of Ar. Because we divide the in- 
tegral by Ar, the will then diverge as Ar approaches 
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Fig. 16. Estimated drift coefficients I)' 1 ' (/i r , r) of the Fokker- 
Planck equation for the road surfaces shown in fig. Q Scales r 
are 108mm (Road 1), 114mm (Road 2), 94mm (Road 3), and 
104 mm (Road 4). Parameterizations are shown as lines. 
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-3-2-1 1 2 3 

Fig. 17. Estimated diffusion coefficients (h r , r) (circles) of 
the Fokker-Planck equation for road surfaces shown in fig. Q 
Additionally the fourth Kramers- Moyal coefficients D^(h r ,r) 
are shown as squares. Scales r are as in fig. 1161 Parameteriza- 
tions are shown as lines. 
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Fig. 18. Estimated drift (a) and diffusion (b) coefficient of 
the Au surface for r = 169 nm. Estimates of D ^ are added as 
squares in (b). 
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Fig. 19. Estimated Kramers-Moyal coefficients for surface 
Crack for a length scale r = 49 /im. (a) D^(h r ,r), (b) 
D (2) (/i r ,r) and D {i) (h r ,r). 



zero. Note that within the same mathematical framework 
the presence of uncorrelated noise can be quantitatively 
determined [33]. 




40 60 
Ar [mm] 



100 




40 60 
Ar [mm] 



100 



Fig. 20. Uncorrelated noise in the case of surface Road 5. 
Dependence of (a) M {1) (h r ,r, Ar) and (b) M {2) (h r ,r, Ar) on 
Ar for r = 104 mm and h r = — 0.60-00. The Markov length 
Im is marked with a dashed line. For illustration a function 
proportional to 1/Ar is fitted to the M (fe) . 



6.2 Conclusions on the estimation of drift and 
diffusion coefficients 

Estimations of the drift and diffusion coefficients (h r , r) 
and D^ 2 \h r ,r) have been performed for all the surfaces 
introduced in section|3 An exception is Road 5, where the 
stochastic process in the scale variable, while still Marko- 
vian, appears to be dominated by additional uncorrelated 
noise. From eqs. (0) and JHJ it can be seen that this leads to 
diverging Kramers-Moyal coefficients I)W , as is the case 
for Road 5. 

As mentioned above, the magnitude of the fourth Kra- 
mers-Moyal coefficient is of particular importance. 
If Z?( 4 ) can be taken as zero, the whole scale dependent 
complexity can be described by a Fokker-Planck equation. 
Otherwise, if is not zero, an infinite set of > is nec- 
essary. In terms of a Langevin equation (|11|) . for ^ 
no Gaussian noise is present. This case is related to un- 
steady stochastic processes [39]. As we see from the to- 
pographies in figs, n and [51 jumps are more likely to be 
present for the Road 4 and Crack surfaces than for the 
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remaining ones. This impression is consistent with the re- 
sult that here we find ^ 0. As a consequence, in these 
cases the Fokker-Planck equation with a drift and diffu- 
sion coefficient is not sufficient to describe the stochastic 
process in the scale variable, because the higher coeffi- 
cients cannot be neglected. The reconstruction of condi- 
tional probabilities (cf. section^) failed for these surfaces. 

The range of scales where the drift and diffusion coeffi- 
cients could be estimated varies for the different surfaces, 
depending on the Markov length on one side and on the 
length of the measured profiles on the other side. In the 
case of Road 2 an additional upper limit for the Markov 
properties was caused by the influence of a strong period- 
icity of the pavement. 



7 Verification of the estimated Fokker-Planck 
equations 

In the previous section methods to estimate the Kramers- 
Moyal coefficients were discussed. We found that this esti- 
mation is not trivial. To prove the quality of the estimated 
we now want to verify the corresponding Fokker- 
Planck equations. 



7.1 Parameterization of Drift and Diffusion 
Coefficients 

With the estimations of the drift and diffusion coefficient 
from section for each surface a Fokker-Planck equation 
jnj is defined which should describe the corresponding pro- 
cess. For the verification of these coefficients it is addition- 
ally desirable to generate parameterizations which define 
D^ x \h r ,r) and D^ 2 '(h r ,r) not only at discrete values but 
at arbitrary points in the (h r , r)-plane. 

Such parameterizations have already been shown in 
figs. El El El and El as lines together with the esti- 
mated discrete values. For it can be seen that for all 
surfaces a straight line with negative slope was used, with 
additional cubic terms for Road 2, Road 4, and Crack. 
The diffusion coefficients were in all cases parameterized 
as parabolic functions. The special shape of the diffusion 
coefficient for Road 2 was parameterized as one inner and 
one outer parabola for small and larger values of h r , re- 
spectively (compare with fig. H7)l . We would like to note 
that both the drift and diffusion coefficients of the cobble- 
stone road presented in [12] are best fitted by piecewise 
linear functions with steeper slopes for larger h r . 

It is easy to verify that with a linear and a con- 
stant the Fokker-Planck equation 10 describes a Gaus- 
sian process, while with a parabolic D^> the distributions 
become non-Gaussian, also called intermittent or heavy 
tailed. For the Au surface it can be seen in fig.Elthat D^> 
has only a weak quadratic dependence on h r and possi- 
bly could also be interpreted as constant (we nevertheless 
kept the small quadratic term because it is confirmed by 
the verification procedure below). If is constant in h r 
the type of noise in the corresponding Langevin equation 



(|ll|l is no longer multiplicative but additive, which results 
in Gaussian noise in the process. Thus the statistics of h r 
in r will always stay Gaussian, and all moments (ft.™) with 
n > 2 can be expressed by the first and second one. As 
a further consequence, the scaling exponents £„ (see sec- 
tion rj are obtained by (h 7 r l ) ~ {h 2 r ) n / 2 as £ n = f£ 2 . This 
linear dependence on n denotes self-affinity rather than 
multi-affinity and is confirmed by the scaling analysis in 
section EHl"! 



7.2 Reconstruction of empirical pdf 

Next, we want to actually evaluate the precision of our 
results. Therefore we return to eq. ©. Knowing and 
£>( 2 ) it should be possible to calculate the pdf of h r with 
the corresponding Fokker-Planck equation. Equation |(SJ| 
can be integrated over ho and is then valid also for the 
unconditional pdf: 



(23) 



(h r ,r) > p{K,r) 



Now at the largest scale Tq where the drift and diffusion 
coefficients could be worked out the empirical pdf is pa- 
rameterized and used as the initial condition for a nu- 
merical solution of eq. J23J. F° r several values of r the 
reconstructed pdf is compared to the respective empirical 
pdf, as shown below in this section. If our Fokker-Planck 
equation successfully reproduces these single scale pdf, the 
structure functions (/i™) can also easily be obtained. 

A second verification is the reconstruction of the condi- 
tional pdf by a numerical solution of Fokker-Planck equa- 
tion © for the conditional pdf. Reconstructing the condi- 
tional pdf this way is much more sensitive to deviations in 
D 1 - 1 ' and D^ 2 ' . This becomes evident by the fact that the 
conditional pdf (and not the unconditional pdf of figs. 1211 
and El determine and and thus the stochastic 
process, see eqs. and ©■ The knowledge of the condi- 
tional pdf also gives access to the complete n-scalc joint 
pdf (eq.^J. Here again the difference from the multiscaling 
analysis becomes clear, which analyses higher moments 
(/i™) = f h™ -p{h r ) dh r of h r , and does not depend on the 
conditional pdf. It is easy to show that there are many dif- 
ferent stochastic processes which lead to the same single 
scale pdf p(h r ). 

For both verification procedures we use a technique 
which is mentioned in [22] and has already been used in 
[12, 14]. An approximative solution of the Fokker-Planck 
eq. © for infinitesimally small steps Ar over which 
can be taken as constant in r, is known [22] 



p(hi,r- Ar\h ,r) 



1 



x exp 



2^7rD( 2 )(/i r , r)Ar 
(tH-hy-D^jh^^Ar) 2 ^ 



(24) 
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A necessary condition for a Markov process is the valid- 
ity of the Chapman-Kolmogorov equation Q21|l [22], which 
allows to combine two conditional pdf with adjacent in- 
tervals in r into one conditional pdf spanning the sum of 
both intervals. By an iterative application of these two 
relations we are able to obtain conditional probabilities 
p{hi, tq — iAr\ho, ro) spanning large intervals in the scale 
r, given that for all involved scales the drift and diffu- 
sion coefficients are known. 

In the following the results of this verification proce- 
dure are shown for those surfaces where the drift and dif- 
fusion coefficients of the Fokker-Planck equation could be 
obtained. 



7.3 Verification results 

The results of the reconstruction of the unconditional pdf 
for the road surfaces with scaling properties are presented 
in fig. The pdf of Road 1 show at smaller scales a 
peak around 5 which is not reproduced by our Fokker- 
Planck equation because in this regime of h r £> (1) and 
-D^ 2 ) could not be estimated with sufficient precision. Here 
it has to be noted that according to eq. (|12J) Coo > ay for 
any r, and thus 5 CToo is a large value for a pdf, denoting 
quite rare events (the r-dependence of ay has been pre- 
sented by S 2 (r) — of, see section 0J). The magnitudes of 
the estimated drift and diffusion coefficients had to be ad- 
justed by a factor of 0.65 to give optimal results in the 
reconstruction. For Road 2 it is likely that the correspon- 
dence between the emprical and reconstructed pdf could 
be improved by a more advanced parameterization of the 
nontrivial shape especially of the estimated drift coeffi- 
cient (see fig. . Here, the estimated drift and diffusion 
coefficients could be used without adjustment. The recon- 
structed pdf for Road 3 are in perfect agreement with the 
empirical ones. A substantial adjustment factor of 0.20 for 
£)W and 0.26 for was necessary to achieve the best 
result. 

Reconstructed conditional pdf are shown in fig. 1221 for 
the road surfaces. While there are deviations for larger val- 
ues of ho, hi, the overall agreement between the empirical 
and reconstructed pdf is good. Especially the rather com- 
plicated shape of the conditional pdf of Road 2 appears to 
be well modelled by our coefficients D^,D^\ As men- 
tioned above, an improved parameterization of may 
lead to even better results. The mag nitudes of L>W,Z)( 2 ) 
were adjusted by the same factors as for the unconditional 
pdf above. 

In the case of the Au film the drift and diffusion coef- 
ficients could be worked out from 281 down to 56 nm, see 
section|Sl In contrast to this regime, the range of correla- 
tion between scales is only about 40 nm, i.e., height incre- 
ments on scales which are separated by at least 40 nm are 
uncorrelated. Nevertheless, both verification procedures 
outlined in section [7] gave good results over the whole 
range from 281 to 56 nm as shown in fig. ESI Here the 
estimated Z)W and D^(h r ,r) were multiplied by factors 
1.3 and 2.2, respectively. 
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Fig. 21. Numerical solution of the Fokker-Planck equation 
12311 compared to the empirical pdf (symbols) for road surfaces 
with scaling properties. For each surface, the topmost solid line 
corresponds to an empirical pdf parameterized at the largest 
scale, and the dashed lines to the reconstructed pdf. Scales 
are (from top to bottom) for Road 1: 316, 158, 79, 66 mm, for 
Road 2: 158, 79, 47, 20 mm, for Road 3: 188, 95, 47, 24 mm. Pdf 
are shifted in the vertical direction for clarity of presentation. 



7.4 Discussion of the verification procedure 

For the verification of the drift and diffusion coefficients 
estimated in section [5] numerical solutions of the Fokker- 
Planck equations © and (|23|l have been performed using 
these estimations. The reconstructed pdf have been com- 
pared to the empirical ones to validate the descripition of 
the data sets as realizations of stochastic processes obey- 
ing the corresponding Fokker-Planck equation. 
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Fig. 22. Numerical solution of the Fokker-Planck equation © 
compared to the empirical pdf for road surfaces with scaling 
properties. Similar to fig. 1121 in each case a contour plot of 
empirical (solid lines) and reconstructed pdf (broken lines) is 
shown on top, with contour levels as in fig. 1121 Below two cuts 
at ho ~ iiJoo are located. Here, empirical pdf are plotted as 
symbols. Scales are ro = 304mm, ri = 158mm (Road 1), ro = 
158 mm, ri = 112 mm (Road 2), and ro = 188mm, n = 92mm 
(Road 3). 
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Fig. 23. Numerical solutions of the Fokker-Planck equations 
Q and compared to the empirical pdf for the Au surface. 

(a) Results of the integrated equation 123H presented as in 
fig. f2"Tl Scales r are 281, 246, 148, and 56 nm (from top to 
bottom). 

(b) , (c) Numerical solution of equation © for the conditional 
pdf compared to the empirical pdf at scales ro = 183 nm, 
ri = 155 nm (b) and ro = 281 nm, n = 56 nm (c). The or- 
ganisation of the diagram is as in fig. 1221 
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Good results were obtained for most surfaces where the 
drift and diffusion coefficients could be derived. In the case 
of Road 4 and Crack we found that the higher Kramers- 
Moyal coefficients and were significantly differ- 
ent from zero, and the empirical pdf could not be repro- 
duced with a Fokker-Planck equation (which only uses 
and _D (2) ). 

It may be surprising that the correspondence between 
the empirical and reconstructed pdf seems better for the 
conditional rather than for the unconditional pdf in some 
cases (compare figs. 1211 and I22fl . One reason may be that 
in fig. [221it is clear that the empirical pdf are not precisely 
defined for combinations of large ho and hi. The eye con- 
centrates on the central regions of the contour plots where 
the uncertainty of the empirical pdf is reduced, as well 
as deviations due to possible inaccuracies and uncertain- 
ties of our drift and diffusion coefficients. This effect is 
also confirmed by our mathematical framework where all 
steps in the procedure are based on the estimation and 
evaluation of the conditional (not the unconditional) pdf. 

The reconstruction procedure allows also to adjust the 
estimated coefficients in order to improve the above-men- 
tioned description, thus compensating for a number of 
uncertainties in the estimation process. While the func- 
tional form of and found in section HJ for all 
surfaces could be confirmed, in most cases the magni- 
tudes of the estimated values had to be adjusted to give 
satisfactory results in this reconstruction procedure. We 
found this effect also when analysing turbulent velocities 
and financial data. One reason may be the uncertainties 
of the estimation procedure. A second source of devia- 
tions may be that the dependence of M^ k \h r ,r, Ar) on 
Ar is not always purely linear in the extrapolation range 
(see section EJ. Thus fitting a straight line and extrapo- 
lating against Ar = may lead to coefficients and 
£)( 2 ) which still have the correct functional form in h r 
but incorrect magnitudes. As mentioned in section in 
our case no general improvements could be achieved by 
the use of nonlinear (i.e. polynomial) fitting functions or 
higher order terms of the corresponding Taylor expansion. 
It is possible that the range of Ar < I m where no Markov 
properties are given is in most cases large enough that ap- 
proximations for small Ar are inaccurate. We would like 
to note that there are also data sets which did not require 
any adjustment of the estimated coefficients, see Road 2 
and [12]. A last remark concerns the latest results in the 
case of Road 1, see fig. Q21 If the fraction of M W (Ar) 
which is proportional to 1/Ar is substracted before per- 
forming the extrapolation, the resulting are substan- 
tially improved in their magnitudes. This may be a way 
to correct the extrapolation of the D^ k > in cases where 
uncorrelated noise is involved. 

In any case, whether an adjustment of and D^ 1 
was needed or not, for the presented surfaces a Fokker- 
Planck equation was found which reproduces the condi- 
tional pdf. Together with the verification of the Markov 
property (TJJ thus a complete description of the n-scale 
joint pdf is given, which was the aim of our work. 



8 Conclusions 

For the analysis and characterization of surface rough- 
ness we have presented a new approach and applied it 
to different examples of rough surfaces. The objective of 
the method is the estimation of a Fokker-Planck equation 

which describes the statistics of the height increment 
h r (x) in the scale variable r. A complete characteriza- 
tion of the corresponding stochastic process in the sense 
of multiscale conditional probabilities is the result. 

The application to different examples of surface mea- 
surement data showed that this approach cannot serve as 
a universal tool for any surface, as it is also the case for 
other methods like those based on self- and multi-affinity. 
With given conditions, namely the Markov property and a 
vanishing fourth order Kramers-Moyal coefficient (cf. sec- 
tional), a comprehensive characterization of a single sur- 
face is obtained. The features of the scaling analysis are in- 
cluded, and beyond that a deeper insight in the complexity 
of roughness is achieved. As shown in [13] such knowledge 
about a surface allows the numerical generation of surface 
structures which should have the same complexity. This 
may be of high interest for many research fields based on 
numerical modelling. 

The precise estimation of the magnitudes of the drift 
and diffusion coefficients for surface measurement data 
still remains an open problem. While for other applica- 
tions a number of approaches have been developed [14,32, 
36, 37] in any case a verification of the estimated Fokker- 
Planck equation is necessary and may lead to significant 
adjustments, as it is the case for some of our data sets. 

We enjoyed helpful and stimulating discussions with R. Friedrich, 
A. Kouzmitchev and M. Haase. Financial support by the Volks- 
wagen Foundation is kindly acknowledged. 
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